Turbulence in Globally Coupled Maps 



M. G. COSENZA AND A. PARRAVANO 
Centra de Astrofisica Teorica, Facultad de Ciencias, Universidad de Los Andes, 
A. Postal 26 La Hechicera, Merida 5251, Venezuela 
(To appear, Phys. Rev. E, 1996) 



OS 



in 



> 

o 
o 
m 
O 

\D 
G\ 

"a 

i 

o 

O 



x 



The phenomenon of turbulence is investigated in the context of globally coupled maps. The local dynamics is given by the 
Chate-Manneville minimal map previously used in studies of spatiotemporal intermittency in locally coupled map lattices. New 
features arise in the globally coupled system; for instance, the transition to turbulence takes place discontinuously at some 
critical values of the parameters of the system. The critical boundaries between different regimes (laminar, turbulent and fully 
turbulent) of the system are calculated on the parameter space. Windows of turbulence are present on some ranges of the 
coupling parameter. The system also exhibits nontrivial collective behavior. A map for the instantaneous fraction of turbulent 
elements is proposed. This map describes many of the observed properties of the system. 



PACS numbers: 05.45. +b, 02.50.-r 



I. INTRODUCTION. 

The transition to turbulence in confined systems and 
its relation to the routes to chaos has been a subject of 
much interest jjj. A very general scenario for the oc- 
currence of turbulence is spatiotemporal intermittency, 
i.e., a sustained regime characterized by the coexistence 
of coherent-laminar and disordered-chaotic domains in 
space and time j^H . There are numerous studies of this 
phenomenon, some of the most extensive of which have 
been on model dynamical systems such as coupled map 
lattices (CML) |,|J||,|,|. The idea is that the ingre- 
dients of a CML — a discrete space, discrete time sys- 
tem of interacting elements whose states vary continu- 
ously according to specific functions — are sufficient to 
capture much of the phenomenology observed in com- 
plex spatiotemporal processes, in particular some rele- 
vant features of spatiotemporal intermittency and tur- 
bulence. In this respect, a CML can be viewed as a 
simplifying replacement for partial differential equations 
of hydrodynamics j9). The onset of turbulence via spa- 
tiotemporal intermittency has mainly been investigated 
on CMLs with local, diffusive-type interactions on both 
low-dimensional Euclidean arrays and on frac- 

tal lattices ||. These studies have permitted the char- 
acterization of the transition to turbulence as a critical 
phenomenon in one and two dimensions and on several 
fractal dimensions with diverse local connectivities. 

Recently, globally coupled elements have been a fo- 
cus of attention in physics and biology [^p]JT^Jl^ . [l^] , 
including systems such as Josephson junction arrays, 
charge density waves, multimode lasers, neural dynam- 
ics, ecological and evolution models. Globally coupled 
maps |]l5| constitute a useful approach to the study of 
many processes on this kind of systems. Globally cou- 
pled maps can be regarded as an extension of a CML 
with diffusive coupling to infinite dimension, or as a lim- 
iting case of a system with long range interactions. Spa- 



tial concepts lose meaning and only temporal properties 
become relevant in globally coupled maps. These charac- 
teristics allow for simpler, mean-field descriptions of the 
behavior of the system. 

In this paper, we investigate the phenomenon of turbu- 
lence in the context of globally coupled maps. The local 
dynamics we employ is the elementary map of Chate and 
Manneville Q, which has been shown to exhibit several 
properties of the transition to turbulence via spatiotem- 
poral intermittency on locally coupled map lattices pL[| . 
Our system of globally coupled Chate and Manneville 
maps provides a situation to compare the roles that local 
and global interactions play on the occurrence of turbu- 
lence. Some new features arise in this model, for instance, 
the onset of turbulence appears to take place discontin- 
uously at critical values of the parameters of the system. 
The model is presented and its statistical properties are 
numerically explored on its parameter space in Section II. 
In section III, a map for the instantaneous fraction of 
turbulent elements in the system is proposed. This map 
model describes much of the behavior observed in sec- 
tion II. The results are discussed in Section IV. 



II. GLOBALLY COUPLED MINIMAL MAPS FOR 
TURBULENCE. 

We consider the globally coupled map system 

N 

xt+i® = (1 - e)/M<)) + /(**(?')), (!) 

i=i 

where Xt(i) gives the state of the lattice element i (i — 
1, . . . , , N = system size) at a discrete time step t, e is the 
coupling parameter, and / is a map describing the local 
dynamics. The instantaneous mean field of the system is 
defined as 
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h t=„52f(xtm (2) 
J'=l 

While the choice of the local map / in Eq. (]l|) de- 
termines the details of the transition to turbulence, it 
need only possess a few general properties to observe 
spatiotemporal intermittency. The map introduced by 
Chate and Manneville Q] contains such minimal require- 
ments, 

(a) f 5(1-|1-2*|), ifx€[0,l] (3) 
J w | x, if x > 1, w 

with r > 2. This map is chaotic as long as f(x) remains 
in the interval [0,1]; however, when f(x) > 1 the dy- 
namics reaches a fixed point in one iteration. The local 
variable can thus be seen as a continuum of stable "lam- 
inar" fixed points (x > 1) adjacent to a chaotic repeller 
or "turbulent" phase (a; G [0,1]). For a coupled system 
the laminar phase is metastable, since sites in this state 
are stable to infinitesimal perturbations but possibly un- 
stable to finite disturbances induced by the coupling. 

In diffusively coupled map lattices, for a large enough 
coupling, the turbulent phase can propagate through the 
lattice in time producing sustained states of spatiotem- 
poral intermittency |^§]. Here, we investigate the corre- 
sponding phenomenon in the context of globally coupled 
maps [Eq. (|l|)] using the local map / [Eq. (j|)]. This situ- 
ation may be considered as spatiotemporal intermittency 
occurring on a diffusively coupled system in a spatially 
high dimensional lattice, or as spatiotemporal intermit- 
tency in a limiting case of long range interactions. Thus, 
the notions of spatial patterns and correlations charac- 
teristic of spatiotemporal intermittency in CMLs with 
local connections are lost and only temporal properties 
of the system become relevant. Above some threshold in 
parameter space (r, e) , and starting from random initial 
conditions, the globally coupled map system relaxes to 
a statistically stationary chaotic regime where each ele- 
ment displays intermittency between laminar and turbu- 
lent phases. As on locally coupled map lattices, the tran- 
sition to this extended chaotic state can be monitored by 
the mean fraction of turbulent sites (F) , a quantity that 
acts as the order parameter for the system H. We have 
calculated (F) as a function of the coupling parameter e 
for several fixed values of r from a time average of the 
instantaneous turbulent fraction Ft, as 

(F) = -j2Ff (4) 

Typically 10 4 iterations were discarded before taking the 
time average in Eq. (^) and r was fixed at the value 
10 4 . The typical system size used in the calculations was 
N = 10 4 . We have found that increasing the lattice sizes 
does not appreciably affect the results presented in this 



paper. The characteristics of the system also persist for 
small lattice sizes (up to N = 100). As initial conditions, 
we use random cell values equally distributed between the 
turbulent and the laminar ranges. It is to be noticed that 
a minimum number of initially turbulent cells is always 
required to reach a sustained state of turbulence. 

Figures 1(a) and 1(b) present plots of (F) versus e for 
two different fixed values of r. For small coupling, the 
system reaches a uniformly laminar state. The onset of 
turbulence takes place discontinuously at a critical value 
of the coupling parameter, e c . The mean turbulent frac- 
tion (F) vanishes at a larger value of the coupling, giving 
raise to a relaminarization of the system; i.e., a second 
transition from turbulence to a uniform laminar state 
takes place at a value e' c > e c , establishing a well defined 
window of turbulence on the range [e c ,e' c ]. Figure 1(b) 
shows than more than one of such windows of turbulence, 
separated by laminar gaps, can occur in globally coupled 
maps as the coupling is varied. 

Windows of spatiotemporal intermittency on a range 
of the coupling parameter have also been observed in cou- 
pled maps on fractal lattices with large enough connec- 
tivities Q. In general, for CMLs with local interactions, 
the transition from a laminar regime to spatiotemporal 
intermittency (and viceversa, in the case of windows of 
turbulence) behaves in many aspects as a second order 
phase transition, characterized by the scaling relation 
(F) ~ (e — e c ) /3 ', where j3 is a critical exponent 
In contrast, the transition between laminar states and 
turbulence in globally coupled maps appears as a discon- 
tinuous jump in the quantity (F), a feature associated 
to first order phase transitions. The absence of spatial 
relations in globally coupled maps rules out the possibil- 
ity of supporting small spatial domains of turbulent cells 
which would be necessary for a continuous transition to 
turbulence. 

The error bars shown on the mean turbulent fraction 
(F) in figs. 1(a) and 1(b) correspond to the standard 
deviation (square root of the variance) of the time se- 
ries of F t at each value of e. With increasing system 
size N, these fluctuations turn out to decrease only up 
to some size, beyond which they remain constant. This 
phenomenon is associated to a nontrivial collective be- 
havior commonly observed in globally coupled maps and 
which has been called "violation of the law of large num- 
bers" 0| : the variance of the temporal fluctuations of 
the mean field [Eq. (□)] does not scale as A -1 for large A, 
but it saturates at some constant value. In the present 
case, the large amplitudes of the standard deviations ob- 
served in figure 1 reflect collective periodic states of the 
system. For example, figure 2 shows the return map of 
the instantaneous turbulent fraction, F t +\ vs. F t for pa- 
rameter values r = 3 and e = 0.295, after discarding 
the transients. The instantaneous turbulent fraction dis- 
plays an approximate periodic behavior, with period six. 
Other nontrivial collective states can be observed at dif- 
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ferent parameter values of the system. The instantaneous 
turbulent fraction is a simpler statistical description than 
the mean field and, in our case, it already manifests a col- 
lective periodic behavior of the system over long times, 
as shown in figure 2. Calculations of return maps of the 
mean field of the system at parameter values giving tur- 
bulence should reveal more details on the nature of those 
collective behaviors, such as global periodic attractors 
and possible quasiperiodic motions. 

We have calculated numerically the critical values e c , 
e' c corresponding to boundaries of the windows of turbu- 
lence observed in the system as a function of r. There 
exists a maximum value r max ~ 3.29 beyond which no 
turbulent regime is observed. In this case, the coupling 
cannot compensate the escape rate from the interval [0, 1] 
in the local maps. Both parameters r and e were varied 
in 1CP 3 in order to detect the regions of turbulence. Fig- 
ure 3 shows the complex critical boundary for the onset 
of turbulence determined in this way. The structure of 
the turbulent windows in parameter space is revealed in 
fig. 3. The fine structure of this critical boundary is ac- 
tually more complex; there are several narrow windows 
of turbulence and laminar gaps at smaller scales. Some 
of the complexity of the critical boundary between tur- 
bulent and laminar regimes has already been suggested 
by less detailed calculations on diffusively coupled one- 
dimensional lattices [|||]. For parameter values inside 
the critical boundary, the globally coupled map system 
exhibits sustained turbulence (0 < (F) < 1). The region 
where full turbulence takes place ((F) — 1) is also shown 
in Fig. 3. In this fully turbulent region, the globally cou- 
pled system, Eq. (pi), uses only the regime x £ [0, 1] of 
the local map, Eq. (m. Thus, the local dynamics is effec- 
tively a tent map with slope r > 2. The global coupling is 
capable of confining all the elements x t (i) to the interval 
[0, 1], even though the local dynamics is repelling in this 
case. As found by Kaneko Jlj, a system of globally cou- 
pled tent maps with r < 2 presents collective behaviors 
(one- and two-band global attractors) and a turbulent 
phase on different regions of its parameter plane (r, e). 
One may expect that a simple extrapolation of the phase 
diagram in |l7| should carry through this structure to the 
range r > 2, corresponding to the fully turbulent region 
in Fig. 3. The laminar regime of the system ({F) = 0) 
occurs for values (r, e) outside the critical boundary. This 
boundary signals the discontinuos transition between the 
two collective states (laminar and turbulent) in the phase 
diagram of the system. Figure 3 also shows a theoretical 
critical boundary for the onset of turbulence in globally 
coupled maps as well as a predicted boundary for full 
turbulence, both obtained from a model presented in the 
next section. 



III. MAP FOR THE TURBULENT FRACTION. 

The existence of two distinct states (laminar and tur- 
bulent) in the local dynamics and the global coupling 
are simplifying features that permit the construction of 
a map for the instantaneous turbulent fraction in the 
system under consideration. In this case, it is possible 
to estimate the instantaneous exchange of cell values be- 
tween the two states in order to compute the change of 
the F t in one iteration. 

Let ((f)// ) be the fraction of cells that being laminar 
(turbulent) at time t will become turbulent (laminar) at 
time t+1. Then, the change in the instantaneous turbu- 
lent fraction from one iteration to the next is 

AF t = F t+1 -F t = <j>+-^. (5) 

A cell i that is laminar at time t becomes turbulent at 
time t + 1 if 

(1 -e)x t (i) + eh t < 1, (6) 

therefore, </>f is the fraction of cells satisfying the condi- 
tion 1 < Xt(i) < x mi where 



1 - e 

The quantity x m is the maximum value that laminar cells 
may have in order to become turbulent in the next iter- 
ation. Thus, the fraction of the laminar range [l,r/2] 
occupied by cells making the transition to turbulence in 
the next iteration is (x m — f)/(r/2 — I). Notice that 

4>t = !: if x m > r/2, . . 

4>+ = o, if x m < 1. [ 1 

For x m £ (l,r/2), we approximate the function (fx// as 

where the positive exponent k\ is a parameter introduced 
to i take into account the non-homogeneous distribution 
of cells in the interval [1, r/2]. 

Similarly, a cell i with Xt(i) < 1/2 becomes laminar in 
the next iteration if 

(1 - e)rx t (i) + eh t > 1, (10) 

while the same transition for a cell with 1/2 < Xt(i) < 1 
occurs if 

(l-e)(l-a: t (t))r + eh t >l, (11) 

therefore, (/>// corresponds to the fraction of cells satis- 
fying the condition x m /r < xt(i) < 1 — x m /r. Notice 
that 
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For x r 
tion <b 



<j>t = 0, if x m > r/2. (12) 
< r/2, we assume the following form of the func- 



k 2 



(13) 



where the quantity (1 — 2x m /r) is the fraction of the tur- 
bulent range [0, 1] occupied by cells that become laminar 
in the next iteration. The positive exponent ki takes into 
account the non-homogeneous distribution of cells in the 
interval [0, 1]. 

By using the above assumptions, the difference map 
Eq. (|J) can be written as 



AF t = 



-F t (1 - 

(1-**) 
I -F t , 



2x m \ 
r / 



k 2 



72-1 



fel 



-F t (l-^) 



if 
if 
if 



X m < 1 
1 *C Xm, 



Figure 4 shows AF t vs. F t obtained from Eq (|l4|) for 
five different values of e and the fixed value r = 2.3. 
The exponents fci = 0.8 and &2 == 1-2 were chosen in all 
the calculations since they optimize the agreement of our 
model with the results from the globally coupled map 
system, but the characteristics of the map model persist 
for a range of values of the parameters k\ and k,2 . The 
real roots (AF t — 0) of Eq (jlij) give the fixed points (sta- 
tionary values) of the turbulent fraction F t . The laminar 
state Ft = is always a stable fixed point. Additionally, 
two other fixed points, F\ and F2 (F\ < F2), occur in the 
interval e S [ec^cL as shown in Fig. 4. The fixed point 
F\ is unstable and it defines a critical value of the in- 
stantaneous turbulent fraction bellow which the system 
invariably falls to the the stable laminar state F t = 0. For 
values of the instantaneous turbulent fraction above this 
^rj.tigal value Fi, the system evolves towards the upper 



fixed point F2, which is stable. As shown in figure 4 and 
— 'in figure 5, the two non- vanishing fixed points coincide 



(14) 



The case x m < 1 (i.e 



0) describes transient tur- 
bulent regimes evolving towards a final laminar state 
(F t = 0). In the case x m > r/2 (i.e cj)// = and 4>f = 1), 
the system reaches a fully turbulent regime (Ft = 1) in 
one time step. 

To obtain the explicit dependence of the difference map 
Eq. ( |l4| ) on the parameters r and e, we need an expression 
for h t . The instantaneous mean field h t can be approxi- 
mated as the sum of two contributions 



h t = h T (r,e)F t + h L (r,e)(l-F t ). 



(15) 



where hx and hx are the instantaneous mean fields of 
the turbulent and of the laminar cells, respectively. In 
particular, the instantaneous mean field hx of the tur- 
bulent cells fluctuates for fixed values of r and e. When 
hx approaches the value 1, the system tends to fall to 
the laminar state; thus, the maximum values of hx in 
parameter space are the relevant ones for the transition 
to turbulence. From the numerical calculations, one can 
roughly approximate hx and hx as: 



hx 



0.48+ 1.05e, 
l + 0.10(r-2). 



(16) 
(17) 



The expression for hx reflects the fact that the range of 
allowed values in the turbulent phase does not depend on 
r, while an increase of the coupling e enhances the flow of 
laminar cells across the phase boundary x = 1, produc- 
ing an increment of the population of turbulent cells with 
x(i) > 1/2. On the other hand, since the range of allowed 
values in the laminar phase is (r/2 — 1), one could expect 
hx to increase with r. The above approximations are suf- 
ficient for the construction of a mapping of the turbulent 
fraction, Eq (|l4|), that describes many relevant properties 
of the globally coupled Chate-Manneville maps. 



at the critical values of the coupling, e c and e' c , corre- 
sponding to tangent bifurcations of the map Eq. (]lj) . In 
figure 5 the fixed points of the map Eq. ( |l4| ) are plotted 
as a function of e, for fixed r = 2.3. Figure 5 clearly 
shows that at the critical values of the coupling, e c and 
e' c , the transition between laminar (Ft — 0) and turbu- 
lent (Ft — F2) regimes predicted by the map Eq. (|lj) is 
discontinuous, as seen in the numerical calculations. The 
mean turbulent fraction (F) obtained from direct simu- 
lations of the globally coupled map system with r = 2.3 
is also shown as a function of e in fig. 5. The stable fixed 
point F2 obtained from the map model agrees well with 
the direct calculation. 

The values e c and e' c corresponding to each value of 
r define the boundary between laminar and turbulent 
regimes in the parameter space. The curves e c (r) and 
e' c (r) obtained from the model are shown in figure 2. By 
comparing with the critical boundary resulting from di- 
rect calculations, one can see that the model given by 
Eq. (|lj) provides a qualitative description of the transi- 
tion to turbulence in the globally coupled map system. 
For values of r close to 2 and small e, the model repro- 
duces quite well the critical boundary. The map model 
also predicts the existence of a boundary for the fully tur- 
bulent region in the parameter plane. On this boundary, 
x m = r/2 and ht = hx', thus from Eq. <M), we get 
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[1 - e(0.48 + l.Q5e) 



(18) 



This curve is also plotted in fig. 2. However, the sim- 
ple model presented here fails to reproduce the complex 
structure of the critical boundary and it predicts a lower 
value for the maximum value of r for which turbulence 
can be sustained. 
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IV. CONCLUSIONS. 



as in bistable maps. 



Some of the features associated to the occurrence of 
turbulence in globally coupled systems have been ex- 
plored in this paper. We have found that the onset of 
turbulence occurs discontinuously at critical values of the 
parameters, as in first order phase transitions. Previous 
studies of diffusively coupled map lattices with Chate- 
Manneville local dynamics have shown that the transi- 
tion to turbulence in those cases is similar to a second 
order phase transition. 

We have made detailed calculations that reveal the 
complexity of the critical boundaries separating the dif- 
ferent regimes of the system. Multiple windows of tur- 
bulence and relaminarization of the system have been 
observed as the coupling parameter is varied. 

Nontrivial collective behavior arises within the turbu- 
lent region in the parameter space of the globally coupled 
system. In fact, the return map of the instantaneous 
fraction of turbulent cells, which is a simple two-state 
statistical description of the system, reflects this collec- 
tive temporal organization. The appearance of nontrivial 
collective behaviors is map lattices with local or global 
couplings is a subject of current research. Much effort has 
been devoted to establishing the necessary conditions for 
the emergence of this behavior. Periodic collective states 
have been found in globally coupled maps belonging to 
some universality class (tent, quadratic, or circle maps) 
0,0. The observation of collective periodic behavior 
in the present system, where the local dynamics has very 
specific characteristics, suggests that this kind of collec- 
tive behavior should be a rather common phenomenon in 
deterministic systems of coupled chaotic elements. 

The proposed map for the instantaneous turbulent 
fraction explains many of the observed aspects of the 
globally coupled Chate-Manneville maps, such as: the 
existence of the critical values of the coupling, e c and 
e' c , and the critical boundary in parameter space for the 
transition to turbulence; the discontinuos character of 
the transition; the existence of an r max ; the presence 
of a threshold for the initial turbulent fraction in order 
to develop sustained turbulence; and the existence of a 
fully turbulent regime and its boundary in the parameter 
plane. 

The global map cannot describe the complex structure 
of the critical boundary and it predicts a lower value of 
r m „. The theoretical critical boundary becomes less ac- 
curate, with respect to the direct simulations, for large 
values of r and e. However, map models for global quanti- 
ties, such as the mean field or other statistical properties 
of the system, constitute a useful approach to the un- 
derstanding of the collective dynamics of coupled map 
lattices. In particular, global maps analogous to Eq. Jl4| ) 
could model the behavior of other globally coupled sys- 
tems whose local dynamics contains two distinct states, 
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FIG. 1. Mean turbulent fraction (F) as a function of the 
coupling parameter e. The local parameter is fixed at (a) 
r — 3; (b) r — 2.6. The error bars indicate ±1 standard 
deviations. 
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FIG. 2. Return map F t +\ vs F t of the instantaneous tur- 
bulent fraction obtained from direct simulation of the glob- 
ally coupled Chate-Manneville maps. Parameters are r = 3, 
e = 0.295; and system size is N = 10 5 . The number of itera- 
tions shown is 10 4 . 

FIG. 3. Phase diagram of the globally coupled map system. 
The different stationary regimes of the system are indicated 
on the parameter plane. The numerically determined critical 
boundaries between laminar and turbulent states and between 
turbulent and fully turbulent regimes are shown with thick 
line. The smooth (thin) curves correspond to the theoretical 
prediction of those same boundaries; upper curve: theoretical 
laminar-turbulent boundary, lower curve: theoretical turbu- 
lent-fully turbulent boundary (Section III). 

FIG. 4. The map AF t vs. F t obtained from Eq (14) plotted 
for five labeled values of e and fixed value r = 2.3. The 
critical values for the onset of turbulence are e c = 0.0145 and 
e' c = 0.488. 

FIG. 5. Numerically calculated mean turbulent fraction 
(F) as a function of the coupling e for fixed r — 2.3. The 
system becomes fully turbulent ((F) = 1) between e ~ 0.31 
and e ~ 0.39. The error bars indicate ±1 standard deviations. 
The stable fixed point F2 (thick line) and the unstable fixed 
point Fi (thin line) of the map Eq. (14) are also shown as 
functions of e for r = 2.3. 
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